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ABSTRACT 

A dispersion  relation  for  surface  plasmons  (retardation 
neglected)  is  obtained  in  the  presence  of  spatial  dispersion. 

The  physical  system  considered  is  one  of  vacuum  separated  by  a 
plane  interface  from  a dielectric  medium  occupying  the  half- 
space  x3  > 0,  and  characterized  by  a nonlocal  dielectric  function 
€(lc||uu  1*3X3) , assumed  to  be  symmetric  in  x3  and  x£.  Here  Icq  is 
a two-dimensional  wave  vector  whose  components  are  parallel  to 
the  interface,  and  uj  is  the  frequency  of  the  electromagnetic 
field  in  the  medium.  The  dispersion  relation  has  the  form 
1 + kdX(k|m)  |x3  - 0,  x3  - 0)  — 0,  and  an  explicit  prescription  for 
obtaining  the  function  X(It||U)  {x3x£)  is  presented.  The  use  of  the 
dispersion  relation  is  illustrated  by  applying  it  to  two  examples: 
(1)  a local  dielectric  constant;  and  (2)  the  nonlocal  dielectric 
function  used  previously  by  Maradudin  and  Mills.  In  both  cases 
previously  obtained  results  are  recovered. 
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Introduction 


Several  years  ago  Ritchie  and  Marusak^  published  a deri- 
vation of  a dispersion  relation  for  surface  plasmons^  in  the 
presence  of  spatial  dispersion.  Their  result  has  the  simple 
form 


<*3 

k2€(k,uu) 


(1.1) 


Here  €(£,u>)  is  the  frequency  and  wave  vector  dependent  dielectric 
constant  of  the  material  out  of  which  the  semi-infinite  medium 
is  composed,  along  whose  planar  interface  with  vacuum  the  surface 
plasmon  propagates.  In  Eq.  (1.1)  kn  is  the  magnitude  of  the 
projection  of  the  three-dimensional  wave  vector  & on  the  plane 
of  the  surface,  and  kg  is  the  component  of  ic  normal  to  the  surface 

A central  assumption  in  the  derivation  of  Eq.  (1.1)  is  that 
the  electrons  in  the  semi-infinite  dielectric  medium  are  reflected 
specularly  from  the  surface.  It  is  this  assumption  which  ultim- 
ately leads  to  the  appearance  of  the  bulk  dielectric  constant 
€(ic,(u)  in  Eq.  (1.1). 

However,  in  any  real  solid  it  is  unlikely  that  all  of  the 
electrons  are  reflected  specularly  from  the  surface.  It  is  more 
likely  that  some  fraction  are  scattered  specularly  and  the  rest 
diffusely,  or  that  more  general  boundary  conditions  obtain.  In 
any  case,  it  would  seem  to  be  desirable  to  have  a dispersion 
relation  for  surface  plasmons  which  does  not  depend,  in  its 
derivation,  on  some  particular  assumption  about  the  nature  of  the 
interaction  of  an  electron  with  the  boundary.  Moreover,  it 
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would  also  seem  to  be  desirable  to  have  such  a dispersion  rela- 
tion which  is  not  based  on  the  assumption  that  the  dielectric 
properties  of  the  medium  supporting  the  surface  plasmon  are  due 
to  a specific  collective  excitation  in  the  medium,  viz.  an 
electron  plasma,  so  that  it  applies  to  systems  in  which 
these  properties  are  associated  with  other  electric  dipole 
excitations  such  as  IR  active  phonons  and  excitons,  for  example, 
or  combinations  of  them. 

It  seemed  to  be  worthwhile,  therefore,  to  try  to  obtain  a 
dispersion  relation  for  surface  plasmons  in  the  presence  of 
spatial  dispersion  of  a form  similar  to  that  of  Eq.  (1.1),  but 
without  invoking  the  restrictive  assumption  of  specular  reflec- 
tion of  electrons  at  the  boundary  of  the  dielectric,  in  which 
the  central  role  is  played  by  the  nonlocal  dielectric  constant 
of  the  medium,  about  which  a minimum  number  of  assumptions  are 
made. 

In  this  paper  we  present  the  derivation  of  such  a dispersion 
relation.  It  is  obtained  in  Section  II,  and  its  use  illustrated 
by  application  to  systems  for  which  the  dispersion  relation  is 
already  known  in  Section  III.  A discussion  of  the  results 
obtained  is  given  in  Section  IV. 


We  assume  a dielectric  medium  which  occupies  the  semi-infinite 
region  xg  > 0.  The  region  x3  < 0 is  occupied  by  vacuum.  Because 
the  system  possesses  infinitesimal  translational  invariance  in 
directions  parallel  to  the  surface,  the  macroscopic  electric 
field  and  the  displacement  in  it  can  be  expressed  in  the  forms 


E(*,t)  - EOWxOe1^  •*■-*«* 


(2.1a) 


D(i,t)  - tf(inu>|x3)ei^ll  * x«  “ 1<Mt  , (2.1b) 

where  2n  - x^+j^kg  , x„  - + *2X2  , and  ^ and  x2  are  two 

mutually  perpendicular  unit  vectors  in  the  plane  of  the  dielectric- 
vacuum  interface.  Within  the  dielectric  medium  the  relation 
between  3(&||U)|x3)  and  fi(l?||U»  |x3)  is  assumed  to  be 

CD 

fl(£|,a»|x3)  - J dxj€(<|«|*3*5)i(<l»|xj)  x3  i 0 , 


x„  as  0 , 


(2.2) 


where  €(£««>  1*3*3)  is  assumed  to  be  symmetric  in  the  variables  x,  and 

— ■ ■■  ■ — 3 


Thus,  we  assume  the  nonlocal  dielectric  tensor  of  the  semi- 
infinite  medium  to  be  isotropic.  This  assumption  is  similar  to 

i 

and,  because  of  the  lower  symmetry  of  the  semi -infinite  medium, 
may  be  as  restrictive  as  the  combined  assumption  by  Ritchie  and 
Marasuk  of  an  isotropic,  nonlocal,  bulk  dielectric  tensor  and 
specular  reflection  of  electrons  from  the  boundary.  The  assumption 
of  isotropy,  however,  can  be  removed  from  both  the  present  theory 
and  that  of  Ritchie  and  Marusak.  The  result  in  each  case  is  a 


ib 


4 O 


\ 


hi  x 


more  complicated  dispersion  relation.  The  assumption  of  specu- 
lar reflection  from  the  boundary,  on  the  other  hand,  is  essential 
in  the  Bitchie-Marusak  approach,  while  it  is  not  in  the  present 
treatment . 


The  Maxwell  equation 


V • 3(x,t)  - 0 

is  now  combined  with  Eq.  (2.2)  and  with  the  equation  of  the 
electrostatic  approximation 

£(£,t)  - - 7cp(x,t)  , 


(2.3) 


(2.4) 


where  <p(x,t)  is  a scalar  potential,  to  yield  an  integro- 
differential  equation  for  (pAwjxg),  the  Fourier  coefficient  of 
<p(x,t): 

m 00 

*1  I ^^(SnujlXgX^^^ijwIx^)  --glj  J dxjeUnmlxgX^)  -gfv  cp(^n»lx3>  “° 


(2.5) 


To  solve  Eq.  (2.5)  we  begin  by  formally  expanding  both 

i 

€(2|(»|xgXg)  and  (d/dXg)6(icBw IX3X3)  in  double  series  Of 

functions  complete  and  orthonormal  in  the  semi— infinite  Interval 

**  % % 

(0,a$ , and  localized  in  the  vicinity  of  the  interface  at  Xg  - 0: 


€<*|H*3*3>  " S 

m,  n-0 


a^tfnuOqp^Xg^x') 


TgJ  €^flu,lx3xP  * ^ ^0bam^Hu,)«Pm(x3)cpn(x3) 


(2.6a) 


(2.6b) 


Although  our  final  result  will  be  independent  of  any  particular 


choice  for  the 


{Vx3>}> 


a convenient  choice,  for  definiteness. 


is  the  set  of  Laguerre  functions  defined  by 

_ (x.)  . 9}.'4Pl3  V*»a> 

vn'x3'  " e nl  » 


(2.7) 


where  ^(*3)  is  the  n^=  Laguerre  polynomial,  and  3 is  a real, 
positive  parameter  with  the  dimensions  of  an  inverse  length. 

Because  of  our  assumption  that  €(k||<i)  jxgx£)  is  symmetric  in 
Xg  and  Xg,  it  is  the  case  that  the  coefficients  |a|Wn( k rtpj )}  are 
symmetric  in  the  subscripts  m and  n. 

We  next  expand  <p(ls||(i)  |Xg)  and  (d/dx^)cp(^uw  |xg)  in  terms  of 
the  {qpm(x3)}: 


so  that 


j dX3<Dn(X3)cp(k||tt)  jXg) 
o 


Bn(^n«»)  - J d*g9a<*3)  dir  ‘ptfiHV 

o 3 


(2.9a 


(2.9fc 


When  we  substitute  Eqs.  (2.6)  and  (2.8)  Into  Eq.  (2.5)  and  use 
the  orthonormality  of  the  {5pn(*3)},  we  find  the  latter  takes  the 


form 


n«0  n«Q  *****  ** 


(2.10 


The  coefficients  (A^jcucp) } and  {Bn(S|joj)3  are  not  independent 
however.  If  we  integrate  by  parts  the  integral  in  Eq.  (2.9b)  we 
obtain  the  result  that 

OD 

Bn(fiw<j))  - - <pn(0)<p(£n<«  10+)  - J dXgcpCtfuu,  lx3)  ^ir  cpn(x3)  * (2-r 

o 

where  the  notation  qp(£||(u|0+)  is  intended  to  emphasize  that  it  is 
the  potential  on  the  medium  side  of  the  interface  that  is  being 
considered.  If  we  substitute  Eq.  (2.8a)  into  the  second  term  on 
the  right  hand  side  of  Eq.  (2.11),  we  obtain  finally  the  relation 


-*n(°)<P(£||«»l0+)  - S em>Vi»tu)  * 


(2.12 


9. 


1 


where 


np 


I “^(cEI  «Pn(^3))<Pp(*3) 

o 3 


(2.13a) 


We  note  for  future  reference  that  e can  be  expressed  equivalently 

np 

as 

CO  CO 

enp  = J to3  J dx3|Pp<x3)  ST  6<x3-x3)cPn(x3>  " " cpp(0)cpn(0)  - epn  . 

o o 3 

(2.13b) 

On  combining  Eqs.  (2.10)  and  (2.12)  we  find  that  the  coeffi- 
cient Bn(fi||uj)  is  given  by 


OB  .J 

Bn(S||0u)  - - kfep(ic1,<M|0+)  ^[kfl  + ea“1B]npcPp(0)  , 


(2.14) 


where  for  convenience  we  have  suppressed  the  argument  (k||(u)  of 
the  matrices  a(£||uu)  and  £(ft||to)  on  the  right  hand  side  of  this 
equation.  It  follows  from  this  result  and  Eqs.  (2.2),  (2.4), 

(2.6a)  and  (2.9b)  that 

D3(£||U)|0+)  " kfcp(2||u>| 0+)  E $ (0)[k^a  1 + ea  1Ba 

m,  n-0 

(2.15) 

We  now  turn  to  the  vacuum  region  x^cO.  In  this  region  the 
function  qp(ic||0u  |x3)  is  readily  found  to  have  the  form 

k||X3 


cp(£|,U)|x3)  - Ae 


x,  < 0 , 


(2.16) 


where  A is  an  arbitrary  constant,  while  the  displacement  compo- 
nent D3(£||U)|x3)  becomes 


d3<2h«m|x3>  * " Ak||® 


k„x3 


x3  < 0 


(2.17) 


10. 


The  boundary  conditions  in  the  problem  are  the  continuity 
of  op(x,t)  and  Dg(£,t)  across  the  surface  Xg  - 0.  The  first  of 
these  yields  the  relation  <p(lc|nu|0+)  - A;  the  use  of  this  result 
together  with  Eqs.  (2.15)  and  (2.17)  in  the  second  boundary 
condition  yields  the  dispersion  relation  for  surface  plasmons 
in  the  form 

m 

1 + kll  E Vm( 0)^(3, jU>)cpn(0)  - 0 , (2.18a) 

m,n-0 

where 

Xm(^«|0>)  - [kfr1  + eT1^1^  . (2.18b) 

It  is  useful  to  point  out  that  is  symmetric  in  m and 

n.  For  this  purpose  it  suffices  to  show  that  the  matrix 

kfa"*1  + ea"  -xb r1  is  a symmetric  matrix.  Inasmuch  as  the  matrix 
*•-1  ** 

a is  symmetric  because  the  matrix  a is,  it  is  necessary  to  show 
only  that  the  matrix  X = esT ^a  1 is  symmetric. 

To  show  this  we  first  relate  the  coefficients  {bmn(^||<»)} 
to  the  coefficients  {amn(£||U>)}  • We  see  from  Eq.  (2.6b)  that 

CD  CO  ' 

ban(*llu,)  ~ I **3  ! dll  6<^llUJ!x3x3><Pn(x3)  * (2-19) 

O O O 

, An  integration  by  parts  yields 

"oo 

bmn(^«<u)  " " J dx3cPm(o)€(^ll«>lox3){Pn(x3>  “ 

o 

00  09 

“ J **3  J <**3  (dir  *m(V)6(*lHx3x3)<*’n(x3')  * (2-20) 

o o ^ 

Ve  now  substitute  Eq.  (2.6a)  into  the  right  hand  side  of  this 
equation,  and  use  Eq.  (2.13a)  to  obtain 


11. 


bnn(*l®>  " “ S {qpm(o)cpr(o)  + emr}  a (£„«») 
r-o 

" S ena  arn(^»tu)  * (2.21) 

r 

The  second  equality  in  Eq.  (2.21)  follows  from  the  use  of  Eq. 
(2.13b). 

It  follows  from  Eq.  (2.21)  that  the  matrix  X can  be  written 
equivalently  as 

X = ea“1(eTa)a“1  = ea“1eT  = ST  , (2.22) 

where  0^  denotes  the  transpose  of  the  matrix  0.  The  matrix 
X(£||<«)  Is  therefore  symmetric,  and  can  be  written  alternatively  as 

^(Sip)  = [k^r1  + ' -(2.18c) 

Equation  (2.18)  may  have  some  interest  of  its  own  in  connec- 
tion with  purely  numerical  studies  of  surface  plasmons  in  spatially 
dispersive  media.  However,  at  this  point  we  direct  our  efforts 
to  re-expressing  Eq.  (2.18)  in  a form  in  which  the  {cp^x^)}-  no 
longer  appear  explicitly  or  implicitly. 

We  begin  by  introducing  the  function  x(£||U>  |xgX^) , which  is 
defined  by 

CO 

X(*H* 3*3)  - 2 cPm(x3)Xmn(^l|(«)«Pn(x^)  x3»x3  > 0*  (2.23) 

m,n=*o 

In  terms  of  this  function  the  dispersion  relation  (2.18)  becomes 

1 + kB  X(£|l<»l00>  = °*  (2.24) 





To  obtain  the  function  ^(£|| tu jx^Xg)  or  more  accurately, 
the  equations  determining  it,  we  return  to  Eq.  (2.18c),  which  we 


rewrite  as 


2 —1  *1 
*1  ^ amp  *pn  + p^s  empapresr*sn  " 6mn* 


(2.25) 


We  now  multiply  both  sides  of  this  equation  by  cpm ( x3 ) cpn  ( 3Cg ) and 
sum  on  m and  n from  zero  to  infinity.  We  recall  the  orthonormality 


and  completeness  conditions 

6mn  * rdX3«Pm(x3)^n(x3) 


(2.26a) 


2 «Pm<X3)<P»(x3)  " 5(x3"  x3K 
m-o 


(2.26b) 


We  next  note  that 


2 'Pm(x3^anm^||a,^«Pn(x3^  ” e (^HU,IX3X3) 


(2.27) 


m,n-o 


where  e“1(S-ti)  |x„xX)  is  the  inverse  of  €(£p<« |x3xi)  in  the  sense 


that 


J ^*3  ^ (Su<w  |x3xg)  € (^|j(u  I X3X3 ) 

- J dx'€  (g,|0)|x3x')€“1(S,|U)|xJx')  - 6(x3-x') 


(2.28> 


We  also  use  the  relations 


^ ^m  ^ x3 ^ emncPn  ^x3  ^ “ 

m,  n-o 


" k&Z  *l*3-*3)  + 6(x3>6(x3)] 


(2.29a) 


- 3!-  6(*3  - *3). 


(2.29b) 
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The  result  of  all  of  this  Is  to  transform  Eq.  (2.25)  into 
kll  ^llU)lx3x3^X^lltulx3x3^  " 

dxT  J ^3  e (^llw  1X3X3> dxT'  X<^tl«»lx3x3)  ” 

w O O 

00 

- 60c3)J  dxj  e'1(8||l.|0x|)^|-,  x<fill»l*3*3>  - S<x3-zp.  (2.30) 

o 3 

Finally,  we  see  that  the  function  X(lc|,<i,|x3x£)  which  satisfies  the 

integro-differential  equation 
00 

k?  1^3  €’1^ll«,lx3x3)x(ill<"lx3x3)  “ 

00 

“ cEtT  J € (^ntt)lx3xg)  ^~»X^||*«!x3x3)  “ .5(x3-Xg),  (2.31a) 

3 O 3 

together  with  the  boundary  condition 
00 

J d*|  €'1(«,»|0  xp  -|y  x<*u«>  1*3*3)  - 0.  (2.31b) 

O 3 

! 

is  a solution  of  Eq.  (2.30) . As  we  are  seeking  electric  fields 
localized  in  the  vicinity  of  the  surface  x3  - 0,  the  boundary 
condition  which  must  be  imposed  on  x(&||tu |XgXg)  at  infinity  is  that 

lim  X(£||t»»|x3x4)  - 0.  (2.31c) 

x3-*° 

Equation  (2.24),  together  with  the  prescription  for  obtaining 
X(£||iu|x3x3)  given  by  Eqs.  (2.31),  is  the  central  result  of  this 
section:  the  surface  plasmon  dispersion  relation.  In  the  next 
section  we  illustrate  its  use  by  applying  it  to  two  cases  for  which 
the  surface  plasmon  dispersion  relation  is  known. 
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III.  Examples 

In  this  section  we  solve  Eqs.  (2.31)  for  two  choices  of 
fXgXg)  , one  local  and  the  other  nonlocal,  for  which  the  surface 
plasmon  dispersion  relation  is  already  known,  to  illustrate  the  use 
of  the  dispersion  relation,  Eq.  (2.24). 

A.  A Local  Dielectric  Tensor 

We  first  consider  the  case  in  which  €(£||<d )xgx3)  is  given  by 

€(£||tu  jXgXg)  - €(uj)6(x3- x3'),  (3.1) 

where  €(tw)  is  independent  of  £||  and  of  the  coordinates  x3  and  x3# 

It  follows  that 

€"1(S„(U|x3x^)  = 6(x3-x').  (3.2) 

The  equation  satisfied  by  x(£nu> JXgX')  in  this  case  becomes 
k2  1 d2 

erfj  1X3X3>  - eftjy  X<gll<«|x3xp  - 8(x3-  x')  (3.3a) 

dx« 


together  with  the  boundary  condition 


X(^„u)|x3x')  I -0 
x3“° 
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(3.3b) 


with  the  aid  of  Eq.  (3.3b),  with  the  result  that 


“kBlx3“x3l 


+ e 


-M*3+*3> 


}• 


(3.5) 


When  Eq.  (3.5)  is  substituted  Into  Eq.  (2.24),  we  obtain  as  the 
surface  plasmon  dispersion  relation 

1 + €(<«)  - 0, 

a well-known  result  for  this  case. ^ 


(3.6) 


r\ 

B.  A Nonlocal  Dielectric  Constant 


The  second  example  we  consider  is  based  on  the  nonlocal 
dielectric  function 

€(t|U> 

I 

which  was  used  by  Birman  and  Sein^  , Maradudin  and  Mills and 
Agarwal,  Pattanayak.,  and.  Wolf^  , in.  studies  of  various  optical  opti- 

cal  properties  ef  semi-infinite  media  in  the  presence  of  spatial 

dispersion.  There  are  some  unphysical  features  of  the 

model  underlying  this  form  for  edc|,aj|x3x^) , as  will  be  discussed 
in  the  following  section.  However,  because  results  based  on  its  use  are 
•▼•liable  it  serves  as  a useful  example  to  illustrate  Eq.  (2.24), 
which  is  not  restricted  by  these  unphysical  features.  In  Eq.  (3.7) 

€Q  la  the  optical  frequency  dielectric  constant,  which  is  assumed 
to  be  real  and  frequency  independent;  is  a plasma  frequency; 

D is  a positive  constant  which  defines  the  curvature  of  the  exciton 
or  transverse  optical  phonon  dispersion  curve  at  & « o;  and  r is 
given  by 


•!  ir|x3-xi| 

" €06(x3-*$)  + 1 e *3»x3  * 0 » <3-7> 


■ 


€ (4Mu)  Ixg*^  “ T~  6(*3”x3)  + 


La 


+ Ae 


The  coefficient  A is  obtained  by  requiring  that  the  expression 
(3.13)  also  satisfy  the  original  integral  equation  (3.9).  In  this 
way  we  obtain  finally  the  result  that 

€ 1(S|ju)|x3x^)  - 6(x3-x')  + 


^ ou 

s 6^x3"*x3^  + f(^||tu  |X3X3)  . 


We  note  the  result  that 
( d2  2^,.* 


(“£Z  + «2)f(gll“>lx3x3)  “ 6(xg-xP 


The  integral  equation  (2.31a)  for  the  function  X(£||W  1*3X3) 


now  takes  the  form 


2*  j ffs 

^ (^2  * *?)*<*•“  1*3*3>  + ^ L *t31<8l“',3x3>  X 


3 “ o 


X x<£|»l*3x3>  - *?  J dxJf(S„«,|x3x')x(g|,tu|x5x')  - 

<1*3  o 

“ -6(x3-x|). 

2 2 2 

We  now  apply  the  operator  (d  /dx3)  + a to  both  sides  of  this 
equation,  and  use  Eq.  (3.15).  to  obtain 


has  been  used 
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Substitution  of  Eq.  (3.20)  into  Eq.  (2.24)  yields  as  the 
dispersion  relation  for  surface  plasmons 

e(ui)  r ikli+0f  ikn-r,  k„ 

1 + l1  + TEpZ  IFjJ+f]  + 2Ir(€((w)  - €o}  “ °‘ 

This  equation  can  be  rearranged  into  the  form 


(3.22) 


1 + €(<u)  - ik„[€(u))  ^ j 


€((«)-€, 


(3.23) 


Inasmuch  as  the  right  hand  side  of  this  equation  is  explicitly 
proportional  to  kN , in  obtaining  the  frequency  of  the  surface 
plasmon  to  lowest  nonvanishing  order  in  kn  it  suffices  to  set 
kn  - 0 in  the  expression  in  braces.  We  thus  obtain  the  approximate 


dispersion  relation 


1 + €((!))  — ik  ||  ^€(ou)  ~ ^ + 


1 IN.  «<•)-«, 


:]k,-o 


(3.24) 


At  this  point  we  neglect  the  intrinsic  damping  in  the  dielectric 
medium  and  let  y tend  to  zero  through  positive  values.  In  this 


limit  we  find  that 


Itk„  - 0)  - iD"*((i>2  - 

- D~*(uu2  - 

«(k|  - Q)  - iD-*(u,2  - 

- D“*(<«2  - 


2v  $ 

(U  )a 


(II  ) * 


(I)  < (II, 


«n  > <u. 


ui  <<u. 


(II  >ui. 


(3.25) 


(3.26) 
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Ve  next  note  that  the  equation  € (uu)  + 1 • 0 is  the  surface 
plasmon  dispersion  relation  in  the  absence  of  spatial  dispersion. 
Consequently  the  frequency  of  the  surface  plasmon  in  the  absence 
of  spatial  dispersion,  uu_,  is  the  solution  of  the  equation 

S 


I 


I 


6(uj8)  +1-0 


Ve  thus  write  the  solution  of  Eq.  (3.24)  in  the  form 


oj  - u>s  + Au>s  , 


(3.27) 


(3.28) 


where  Auug  is  of  first  order  in  k||  . Substituting  Eq.  (3.28)  into 
Eq.  (3.24)  and  linearizing  the  resulting  equation  in  Auu_  with  the 
aid  of  Eq.  (3.27)  we  obtain 


„ lk"  fl  3*o1 
€#(cO  L®  2T  j 


*11  “0 
uu  - uu* 


(3.29) 


If  we  denote  the  static  (uu-0)  value  of  the  bulk  dielectric 


constant  at  k - 0 by  € . we  have  the  relation 

s 

6s  * €o  + (<"p/%)  * 


(3.30) 


With  the  substitution  of  Eq.  (3.21)  (with  y“0)  into  Eq.  (3.27)  we 
find  that  the  frequency  of  the  surface  plasmon  in  the  absence  of 


spatial  dispersion  is  given  by 


2 2 % 
s "o  T l+€. 


(3.31a) 


or  equivalently 
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It  is  readily  determined  that 

% < “s  < “L  * 


(3.31b) 


(3.32) 


These  inequalities,  together  with  Eqs.  (3.25)  and  (3.26),  enable 

us  to  express  Atu  given  by  Eq.  (3.29)  in  the  form 

s 


°ikn  r i w0+*>  i ~] 

€'(u)s)  L C^-ouf)*  2 (u>2-uj2)4  -I 


D*(l +€  )*k 


V 


go>ZkH  r 
'<u>J  L 


* _ 


-]■ 


(3.33) 


which  agrees  completely  with  the  result  for  obtained  by 

s * 

(5\ 

Maradudin  and  Mills  on  letting  the  speed  of  light  tend  to 
infinity  in  the  surface  polariton  dispersion  relation  they 
obtained  on  the  basis  of  the  dielectric  constant  (3.7). 
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IV.  Discussion 

In  this  paper  we  have  obtained  an  explicit  dispersion  relation 
for  surface  plasmons  at  the  interface  between  vacuum  and  a semi- 
infinite dielectric  medium  occupying  the  half-space  xg > 0 and 
characterized  by  a nonlocal  dielectric  constant  € (Sj|  uu|  XgXg)  . It 
should  be  remarked  that  the  derivation  of  Eqs.  (2.24)  and (2.31) 
presented  here  is  somewhat  artificial.  It  should  be  possible  to 
obtain  these  results  more  directly,  but  we  have  not  attempted  to 
do  so.  It  also  appears  as  if  the  use  of  Eqs.  (2.24)  and  (2.31) 
to  obtain  the  surface  plasmon  dispersion  relation  may  involve 
lengthier  calculations  than  are  required  in  a more  conventional 
approach  based  on  a direct  solution  of  the  equations  of  electro- 
statics. Nevertheless,  it  still  seems  useful  to  have  an  explicit 
dispersion  relation,  Eq.  (2.24),  with  a definite  prescription  for 
obtaining  the  function  x(£(|«) | *3X3)  entering  it,  Eq.  (2.31).  For 
example,  in  the  case  of  dielectric  constants  of  more  complicated 
form  than  those  considered  in  the  preceding  section  Eq.  (2.31) 
can  serve  as  the  basis  for  approximate  determinations  of  X(lt|(u>| XgXg) , 
e.g.  by  variational  methods. 

The  fact  that  the  method  developed  here  yields  the  same  surface 
plasmon  frequency,  Eqs.  (3.28),  (3.31),  and  (3.33),  as  was  obtained 

( 5) 

by  a rather  different  method  by  Maradudin  and  Mills'  7 is  of  some 
independent  interest  for  the  following  reason.  The  dielectric 
constant  (3.7)  is  obtained  by  partially  Fourier  transforming  the 
bulk  dielectric  constant 


€ (k,w) 


2 

■d 

~2 5* — 2 

uu*  + Dk  - oj  - iiuy 


(4.1) 


\ 


) J 

1 

1 

I j 

I 


/ (7 ) 

and  restricting  xg  and  xg  to  be  positive.  Recently  Horingv,/  used 
£(£,<*>)  as  given  by  Eq.  (4.1)  in  the  Ritchie-Marusak  dispersion 
relation,  Eq.  (1.1),  and  obtained  an  expression  lor  the  frequency 
of  a surface  plasmon  which  differs  from  the  Maradudin-Mills  result, 
Eqs.  (3.28),  <3.31),  and  (3.33).  This  difference  is  not  very 
surprising  in  fact.  As  we  have  noted  in  the  Introduction,  under- 
lying the  Ritchie-Marusak  dispersion  relation  is  the  assumption 
that  electrons  are  reflected  specularly  from  the  surface  of  the 
solid.  The  dielectric  constant  (3.7),  on  the  other  hand,  repre- 
sents the  so-called  dielectric  approximation,  in  which  the  bulk 
dielectric  constant  is  partially  Fourier  transformed  according  to 
Eq.  (4.2),  after  which  x^  and  x^  are  restricted  to  be  positive,  so 
that  the  resulting  € (ic||U)  |xgXg)  is  no  longer  a function  of  Xg-Xg, 

but  depends  on  Xg  and  Xg  separately.  It  is  known  that  the  dielectric 

(8) 

approximation  does  not  conserve  particle  number v ' , and  causes  the 

(Q) 

surface  to  act  as  a source  or  of  energy'  ' , neither  of  which 

is  the  case  when  specular  is  assumed.  Thus  a different 

physical  situation  is  being  considered  in  Horing's  work  from  that 
assumed  in  the  work  of  Maradudin  and  Mills.  Consequently,  it  is 
not  surprising  that  different  surface  plasmon  frequencies  are 
obtained  in  these  two  calculations.  However,  this  example  points 
up  the  desirability  of  having  a dispersion  relation  for  surface 
plasmons  which  is  valid  for  more  general  physical  situations  than 


is  represented  by  the  assumption  of  specular  reflection  at  the 
surface.  Equations  (2.24)  and  (2.31)  provide  such  a dispersion 
relation. 
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